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Recent measurements of solar wind turbulence report the presence of intermittent, exponentially 
distributed angular discontinuities in the magnetic field. In this Letter, we study whether such 
discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the dis- 
continuities by measuring the fluctuations of the magnetic field direction, A9, across fixed spatial 
increments Ax in direct numerical simulations of MHD turbulence with an imposed uniform guide 
field Bo- A large region of the probability density function (pdf) for AO is found to follow an expo- 
nential decay, proportional to exp(— A9/0,), with characteristic angle 9, w (14°)(6 rms /_Bo) ' 65 for 
a broad range of guide-field strengths. We find that discontinuities observed in the solar wind can 
be reproduced by MHD turbulence with reasonable ratios of b Yms /Bo- We also observe an excess of 
small angular discontinuities when Ax becomes small, possibly indicating an increasing statistical 
significance of dissipation-scale structures. The structure of the pdf in this case closely resembles the 
two-population pdf seen in the solar wind. We thus propose that strong discontinuities are associ- 
ated with inertial-range MHD turbulence, while weak discontinuities emerge from dissipation-range 
turbulence. In addition, we find that the structure functions of the magnetic field direction exhibit 
anomalous scaling exponents, which indicates the existence of intermittent structures. 



PACS numbers: 95.30.Qd, 96.50.Ci, 96.50.Tf 

Introduction. — Over the past several decades, multiple 
spacecraft have measured the fluctuations of the mag- 
netic and velocity fields in the solar wind [TJ This 
has provided a wealth of data from which to test plasma 
models such as magnetohydrodynamic (MHD) turbu- 
lence [31 0]. Among other things, this allows one to 
study the intermittency of MHD and dissipative struc- 
tures such as current sheets. However, the question re- 
mains of whether the solar wind can be adequately de- 
scribed by MHD turbulence, and whether different types 
of structures can be distinguished from the data. 

The existence of intermittent structures in a plasma is 
implied by abrupt changes in magnetic field directions. 
Bruno et al. [5] studied this feature in solar wind data 
from the Helios 2 spacecraft by performing a minimum 
variance analysis on the magnetic field vector. It was 
found that there are times when the magnetic field un- 
dergoes large changes in direction, implying that the solar 
wind contains intermittent structures. It was proposed 
that these strong fluctuations mark the boundaries of 
flux tubes that originate in the Sun and are passively 
advected by the solar wind. 

Later, Borovsky [6] used data from the ACE space- 
craft to examine the probability density function (pdf) 
of angular shift in the magnetic field, given by A9 = 
cos -1 [B\ ■ B2/I-B1H-B2I) where B\ and B2 are mea- 
surements of the magnetic field taken at two different 
times. Two populations of magnetic discontinuities were 
discerned. The first population consists of strong discon- 
tinuities at 30° < AO < 170° with an exponentially de- 



caying pdf proportional to exp(— A0/24.4 ). The second 
population consists of weak fluctuations at 5° < A9 < 
30° which can be fit by exp(-A0/9.4°). Miao et al. [7 
performed a similar analysis on slow wind data from the 
Ulysses spacecraft. In their case, the first population 
spanned approximately 50° < AO < 160° and was pro- 
portional to exp(— A0/3O.O ), while the second popula- 
tion spanned 30° < Ad < 50° and was proportional to 
exp(— A0/18.6 ), broadly agreeing with Borovsky's re- 
sult. In both reports, the strong discontinuities were in- 
terpreted to come from coronal flux tube walls, while the 
weak discontinuities were assumed to be turbulent fluc- 
tuations. 

An alternative hypothesis is that the magnetic discon- 
tinuities in the solar wind are predominantly generated 
by nonlinear interactions [BTfTU] . In this case, the discon- 
tinuities evolve dynamically as the solar wind expands 
away from the Sun. It is known that current sheets form 
spontaneously in MHD turbulence, providing a natural 
source of discontinuities. 

These studies raise a principle question of why the 
pdfs of angular shifts have exponential laws (regardless 
of their origin), and what determines the typical angu- 
lar discontinuities characterizing the two observed scal- 
ings. Although the statistical properties of magnetic 
fields in MHD turbulence have been studied before [TlT - 
IT3"] , the statistical properties of magnetic discontinuities 
have been addressed to a lesser extent [14j [15] . In par- 
ticular, the statistical properties of angular shifts have 
not been studied in numerical simulations of MHD tur- 
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bulence. In this Letter, we investigate the statistical 
properties of magnetic discontinuities in direct numeri- 
cal simulations of MHD turbulence with an imposed uni- 
form magnetic field B . We use simulations with sev- 
eral choices of So that span the transition from turbu- 
lence with a weak mean field to turbulence with a strong 
mean field. We find that the pdfs of A8 all contain a 
region well fit by P(A9) oc exp(— AO/ 8*), with a scal- 
ing that depends on the background field strength as 
0» w (14°) {b lms /B ) 0M5 . The pdf of strong discontinu- 
ities in Borovsky's and Miao et al.'s studies could then 
be explained by MHD turbulence with reasonable fluc- 
tuations to guide-field ratios for the solar wind: the first 
population in Borovsky's results is consistent with our 
results for b Tms /B ~ 2.4, and Miao et al.'s result is con- 
sistent with b lms /Bo ~ 3.2. 

We also find that when the spatial increment Ax be- 
comes small enough, the pdf exhibits an excess of small 
angular discontinuities, which closely resembles the "sec- 
ond" population in the solar wind observations. We con- 
jecture that this population indicates the increasing con- 
tribution of structures near the dissipative scale, where 
the nature of turbulence changes. 

One can then envision the following explanation of the 
solar wind observations in which both populations of dis- 
continuities arise from MHD turbulence. According to 
our analysis, the strong discontinuities may be produced 
by discontinuities in MHD turbulence with b lms > B , 
which is typical of the solar wind. These discontinuities 
will be observed regardless of whether any flux tubes are 
advected by the solar wind. The population of weak dis- 
continuities, on the other hand, may be produced by a 
different kind of turbulence, which still generates an ap- 
proximately exponential pdf but with a different char- 
acteristic angle. In particular, this population may be 
associated with near-dissipation-scale turbulence. If so, 
this population emerges only when the interval between 
measurements is small enough. 

This picture is attractive because it explains the origin 
and the exponential pdf of large angular discontinuities 
as well as the small ones. This model for the solar wind 
magnetic discontinuities can be complementary or possi- 
bly alternative to Borovsky's flux tube hypothesis. 

Method. — The incompressible MHD equations can be 
written as 

d t v + (v ■ V)v = -Vp + (V x B) x B + vV 2 v + f, 
d t B = V x x B) + r]W 2 B, 
V • v = 0, 

V-B = 0, (1) 

where v(x,t) is the plasma velocity, B(x,t) is the mag- 
netic field, p is the pressure, and f(x,t) is the external 
forcing. We take the viscosity v and resistivity r\ to be 
equal. 
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FIG. 1. The pdfs for angular shifts in magnetic field for MHD 
turbulence with fixed guide field Bo — 1 and measurement in- 
crements in the range 2 < Aa; < 32. Each pdf has a region 
that is approximately an exponential decay, with a charac- 
teristic angle that is independent of Aa:. Note that for Aa; 
near the dissipation scale, small angular shifts become more 
abundant than expected from the exponential tail. 



The angular shift in magnetic field between two space- 
craft measurements is the angle between B(t) and B(t + 
At), where B(t) is the magnetic field vector at time t 
and At is the time increment. Instead of time incre- 
ments, we use a spatial increment in our analysis, which 
can be approximately related to the time increment by 
Ax ps VswAt, where Vsw is the solar wind velocity. 
The angular shift in magnetic field between two points 
P\ = (x, y, z) and P2 = (x + Ax, y, z) is then given by 
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cos 



B(x,y,z) ■ B(x + Ax,y,z) 
\B(x,y,z)\\B{x + Ax,y,z)\ 



(2) 



We analyze data from simulations of driven incom- 
pressible MHD with five different guide-field strengths: 
-Bo G {0.25, 0.5, 1, 5, 10}, in comparison to the root mean 
square average perpendicular fluctuations of b rms ~ 1.3. 
The simulations solve the full MHD equations with a 
Reynolds number Re « 2200 and conditions similar to 
Ref. [T5]. The simulations have a resolution of 512 3 . 

For each case of Bq, multiple snapshots of the simula- 
tion during steady-state are analyzed. In each snapshot, 
we measure A8 on a line of points separated by Ax in 
the x-direction (this is perpendicular to the guide field, 
but the results in the parallel direction are similar). 
Spatial increments are taken in the range 2 < Ax < 32, 
where Ax — 1 is equal to the mesh spacing. This is 
repeated along other lines in the snapshot with a spacing 
of Ay and Az grid points, where Ay and Az are varied 
to give enough statistics. For each choice of Ax and B , 
we measured A8 at approximately 3.5 x 10 6 data points. 

Results. — The pdf for angular shifts in magnetic field 
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FIG. 2. The pdfs for angular shifts in magnetic field with 
fixed Ax = 8 and several background magnetic fields (where 
Ls ~ 1.3). Each pdf has a region where it takes the form of 
exp(— A9/9 t ), where 9, depends on Bq. 



is shown in Fig. [T] for fixed guide field of Bo — 1 and 
spatial separations in the range 2 < Ace < 32. The axes 
of the plot are log-linear, so a straight line represents 
an exponential function. Each pdf has a region that is 
fit by exp(— A0/6**), where the characteristic angle 6* = 
19.4° is independent of Ax. For small spatial separations 
(Ax < 4), the exponential region is limited to the tail of 
the pdf, and small angular shifts are more abundant than 
expected from the exponential fit. This deviation from an 
exponential can be attributed to dissipative effects, and 
provides a natural explanation for a distinct population 
of weak angular discontinuities. 

Next, we consider how the pdf for angular shifts 
changes when the guide field is varied. This is shown 
in Fig. [2] for fixed Ax = 8 and varying B . Each 
case has an exponential region, but with characteris- 
tic angle 6** decreasing with stronger guide fields. The 
value of is acquired by measuring the slope in the 
exponential region with a least-squares fit. The re- 
sult is that B = [0.25, 0.5, 1, 5, 10] correspond to 6* = 
[43.3°, 34.8°, 19.4°, 6.1°, 3.3°], respectively. 

A functional form of the relation between 8* and 
Bo/brms is found by making a fit to our numerical mea- 
surements. The best fit is a power law (Fig. [3]), given 
by 

9, « (14°) ( WBo) - 65 . ( 3 ) 

It is interesting to point out that the mean angular 
shift obeys a different empirical relation (obtained by fit- 
ting the simulation data), 

(AO) = J AdP{Ad)dA6 cx arctan ^ L1 ^ rms ^ , (4) 
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FIG. 3. The characteristic angle measured in the exponen- 
tial part of the pdf of A9, denoted 9*, is found to have a 
power law dependence on Bo/b T ms- The fit shown is 9* ~ 
(14 )(6 rms /_Bo) 0,65 - It is interesting to note that the mean 
angular shift, (A9), instead has an arctan dependence, shown 
in the inset. 

which is shown in the inset of Fig. [3j This reflects the 
fact that a significant contribution to the average comes 
not from the exponentially declining tail of the probabil- 
ity density function, but rather from its bulk part that 
has a different scaling. Hence, the mean value of AO is 
determined by fluctuations outside of the exponential re- 
gion. In support of this view, we note that the location 
of the maximum in the function A9P(A9) indeed has an 
arctan dependence on Bq. This result reflects the fact 
that a strong background magnetic field prevents large 
deviations of the magnetic field vector from B = BqZ. 

In addition to pdfs, another statistical tool used to 
study intermittency is the spatial scaling of structure 
functions [17l [18] . While the anomalous scaling of var- 
ious structure functions has been observed in fluid and 
in hydromagnetic turbulence, the intermittency of mag- 
netic angular discontinuities has not been studied before. 
We define the structure function of order n for magnetic 
field direction to be 

S$(Ax) = (\M{*x)\ n ), (5) 

where the brackets (• • • ) indicate the spatial average of 
the enclosed quantity. We assume that structure func- 
tions have the form 

S%(Ax) cx (Ax) c ", (6) 

where are the exponents of the structure functions. 

We determine („ by measuring the slope of Sg ver- 
sus Ax on a log-log plot using a least-squares fit. As 
shown in Fig. [4j the relationship between £„ and n is 
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FIG. 4. The structure function exponents, £ n , associated with 
the structure functions of magnetic field direction, Sg . The 
deviation from a straight line implies the existence of inter- 
mittent structures . 



not linear, which indicates anomalous scaling due to in- 
termittency. Flattening around Cn ~ 1 implies the pres- 
ence of shocklike magnetic angular discontinuities, that 
is, current sheets [19j|2Qj. It is important to note that 
measurements of structure functions above order 5 or 6 
in a finite series of data may be contaminated by contri- 
butions from rare events |21) , hence the rightmost values 
in Fig. [4] could be unreliable. 

In addition to angular shifts in the magnetic field, one 
may consider studying discontinuities in velocity field. 
However, the angular shift in velocity is not invariant 
under a Galilean transformation, so it has an ambigu- 
ous physical interpretation. Instead, consider the veloc- 
ity jump, 



At; = \v(x + Ax, y, z) — v(x, y, z)\ 



(7) 



which also has a pdf with an exponential region (not 
shown) (T7]- A scatterplot not shown here reveals that 
AO is correlated with Av at large values, as long as Ax 
is chosen to be small (Ax < 8). 

Interestingly, a similar correlation was found in solar 
wind data for fractional velocity jump Av/v instead of 
velocity jump, where v is approximately constant (the 
solar wind velocity) [5J. The correlation between Av/v 
and AO at large values of AO was attributed to the 
crossing of flux tube walls, but our results suggest that 
this can be explained by current sheets produced from 
turbulence. 

Conclusions. — We have addressed to what extent in- 
compressible MHD turbulence can describe the magnetic 
discontinuities in the solar wind. This was done by study- 
ing the statistical properties of angular shifts in turbu- 
lent MHD simulations with varying strengths of guide 



field. We found that the pdf of angular shifts has a re- 
gion of exponential decay, P(A0) cx exp(— AO/6*), in 
agreement with the solar wind. We found that the as- 
sociated characteristic angle is independent of spatial in- 
crement Ax, but depends on the guide-field strength as 
0, « (14°)(fe rms / J B ) - 65 . 

Our results have implications for the discontinuities ob- 
served in the solar wind. The exponential pdf of strong 
discontinuities in the solar wind is consistent with MHD 
turbulence for values of b rms /Bo that are typical of the 
solar wind. Additionally, our simulations show an excess 
of weak discontinuities when spatial increments are small, 
which can be interpreted as a second population. Turbu- 
lence can therefore qualitatively reproduce the solar wind 
observations, with the strong discontinuities associated 
with inertial-scale turbulence, and the weak discontinu- 
ities with near-dissipation-scale turbulence. MHD turbu- 
lence can then be considered a complementary or alterna- 
tive explanation of the discontinuities in the solar wind. 
Further support for this picture comes from the correla- 
tion between large velocity jumps and angular shifts in 
magnetic field observed in our simulations. 

We also found that the structure functions of the mag- 
netic field direction provide a means to study intermit- 
tency. The observations of magnetic discontinuities in 
the solar wind, and their connection to MHD turbulence, 
motivate further study. A more detailed analysis of these 
results will be presented in a future paper. 
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